N A » ~ s + 上 1\ 
内 燃 机 排 气 余热 回收 梯级 分 段 温差 电 单 偶 模 拟 分 析 
机 琦 ， 舒 歌 群 ， 田 华 ， 卫 海 桥 ， 梁 兴 雨 ， 孙 秀 秀 
(天 津 大 学 内 燃 机 燃烧 学 国家 重点 实验 室 ， 天 津 ，300072) 

摘 要 : 本 文选 取 了 BizTes 和 CoSbs 两 种 温差 电 材料 设计 了 一 种 梯级 分 段 温差 电 单 偶 ， 并 对 其 建立 了 数学 模型 ， 导 出 温差 电 
单 偶 的 功率 和 效率 计算 公式 ,分 析 冷 热 端 陶瓷 片 表 面 温度 、 温 差 电 单 偶 长 度 以 及 表面 对 流传 热 系 数 对 分 段 温差 电 单 偶 性 能 的 影 
响 。 分 析 结 果 表明 提高 热 面 温 度 、 降 低 冷 面 温度 和 提高 热 表面 对 流传 热 系数 均 可 以 提高 分 段 温 差 电 单 偶 的 最 大 输出 功率 和 最 大 
转换 效率 。 缩 短 温差 电 单 偶 的 长 度 可 以 增 大 最 大 输出 功率 ， 但 会 导致 最 大 转换 效率 的 明显 下 降 。 
关键 词 : 内 燃 机 ; 排 气 ， 余热 利用 ; 温差 发 电 
中 图 分 类 号 : TK402 ”文献 标志 码 : A 


Simulation and Analysis of Segmented Thermoelectric Unicouple Used for 


Waste Heat Recovery of lnternal Combustion Engine 
JIA Qi SHU Ge-Qun Tian Hua Wel Hai-Qiao Liang Xing-Yu Sun Xiu-Xiu 
(State Key Laboratory of Engines, Tianjin University, Tianjin 300072, China) 
Abstract: In this paper, a segmented thermoelectric unicouple is designed using two kinds of thermoelectric 
materials--BizTes and CoSb;. The mathematical model of the unicouple is established. Based on the mathematical 
model, calculation formula of the output power and conversion efficiency are derived. And the thesis analyzes the 
influence of surface temperature of hot and cold side ceramics, thickness of unicouple and convection heat transfer 
coefficient of hot surface on the performance of segmented thermoelectric unicouple. The results show that increase 
hot surface temperature, reduce cold surface temperature and increase the convection heat transfer coefficient can 
improve the maximum output power and conversion efficiency of the thermoelectric unicouple. Cutting down the 
thickness of unicouple can increase the maximum output power, but it will result in the reducing of maximum 
conversion efficiency. 
Key Word: internal combustion engine; exhaust; waste heat recovery; thermoelectric generator 
人 采用 分 段 结构 ， 将 适用 于 不 同 温度 的 热电 材料 
沿 温度 梯度 串联 ， 使 单 偶 中 每 段 材 料 均 在 最 佳 的 温 
内 燃 机 输出 的 有 效 功 率 一 般 只 占 燃油 燃烧 总 热 ” 度 内 工作 ， 可 以 提高 温差 发 电器 在 大 温差 条 件 下 的 
量 的 30% 一 45%( 柴 油 机 ) 或 20% 一 30% (汽油 机 )， ”转换 效率 四 ,因此 , 为 满足 对 内 燃 机 排 气 余热 的 回收 ， 
而 其 余 大 部 分 的 能 量 都 通过 冷却 水 和 排 气 散 失掉 ”本 文 设计 了 一 种 分 段 温差 电 单 偶 , 高 温 侧 选 用 CoSb3 
了 ， 造成 能 量 的 严重 浪费 门 。 温差 发 电 是 一 种 绿色 环 。 ” 基 温 差 电 材料 ， 低 温 侧 选 用 BizTes 基 温 差 电 材料 ， 
保 技术 ， 有 具有 性 能 稳定 、 无 噪音 、 体 积 小 、 使 用 寿 。 并 对 其 进行 模拟 分 析 ， 计 算 不 同 热 面 温 度 、 冷 面 温 
命 长 等 优点 。 内 燃 机 余热 源 的 分 析 结果 显示 ， 由 了 度 、 温 差 电 单 偶 长 度 以 及 热 表 面 对 流 传 热 系数 对 其 
排 气 温度 较 高 ， 与 环境 温度 间 存 在 很 大 的 温差 ， 在 ”性 能 的 影响 。 
排 气 系统 使 用 温差 发 电 技术 回收 热量 的 潜力 很 大 中。 、 二 
目前 生 中 日 活 心 vs 证 [si 1 计算 模型 
前 限制 温差 发 电 广泛 应 用 的 主要 原因 是 相对 
较 低 的 热电 转换 效率 。 温 差 发 电器 的 效率 主要 取决 如 图 1 为 温差 电 单 偶发 电 示 意图 ， 每 个 温差 电 
于 温差 电 材料 的 优 值 系数 ， 冷 热 端 的 温度 以 及 温差 ”” 单 偶 主要 由 了 型 和 NN 型 电 偶 臂 ， 以 及 冷 热 两 端 用 于 
电 单 偶 的 长 度 和 截面 积 。 温 差 电 材料 的 优 值 系数 ” 导热 绝缘 的 陶瓷 片 组 成 。 本 文 设计 的 温差 电 单 偶 所 
ZT=a?oT/K, 其 中 a、o、k 分 别 为 材料 的 塞 贝克 系数 、 ”选取 的 P 型 高 温 材 料 为 Bao3Ino3FeCo3Sbi2、 低 温 侧 
电导 率 、 热 导 率 31, 材料 的 优 值 系数 越 大 ， 在 相同 的 为 Bio agsSb1 s2Te; N 型 高 温 材 料 为 Bao4Ino4Co4Sb1>、 
温差 下 ， 热 电 转换 效率 也 就 相对 越 高。 一 般 来 说 ， 低温 侧 为 BizTez7Seo3, 材料 物性 参数 参考 文献 说。P、 
任意 一 种 温差 电 材料 只 会 在 一 定 的 温度 范围 内 取得 。” N 型 电 偶 辟 高 温和 低温 材料 的 长 度 比例 51:6: 为 1:1。 
较 高 的 优 值 系数 ， 严 重 影 响 了 温差 发 电 在 内 燃 机 排 ”假定 热量 只 沿 着 电 偶 臂 方 向 传递 ， 忽 略 电 偶 臂 和 环 
气 余热 这 种 大 温差 余热 回收 条 件 下 的 转换 效率 。 境 的 传 热 损 失 ， 忽 略 接 触 热 阻 和 接触 电阻 。 温 差 电 


| 


作者 简介 : 贾 琦 (1989--)， 男 ， 人 硕士 研究 生 ， 主 要 从 事 内 燃 机 余热 


I 


收 研究 。 通 信 作 者 : 田 华 ， 副 教授 ，thtiu@tju.edu.cn 


201802.00111v1 


chinaXiv 


单 偶 的 上 下 两 面 分 别 是 热源 和 冷 源 ， 在 冷 热源 的 共 
同 作用 下 ， 热 端 陶 次 片上 端面 温度 《以 下 简称 热 面 
温度 ) 为 7,， 冷 端 陶 疙 片 下 端面 温度 (以 下 简称 冷 
面 温度 ) 为 7.,， 电 偶 辟 两 端的 温度 分 别 为 Th 和 7.。 


图 1 分 段 温差 电 单 偶发 电 示意 图 


Fig.l Schematic diagram of segmented thermoelectric 


unicouple 


图 2 为 P 型 电 偶 臂 两 材料 分 界面 温度 的 计算 流 
程 ， 首 先 初 设 分 界面 温度 Tip， 计 算得 到 | ”an(rj7 


和 上 X(T7H7 ， 其 中 (站 和 2 人 (站 分 别 为 高 温和 


低温 材料 的 导热 系数 随 温度 变化 的 函数 。 在 稳 态 过 
程 中 ， 高 温和 低温 材料 垂直 于 电 偶 怖 方向 的 截面 上 


he 
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导热 量 相等 ， 故 截面 上 的 导热 量 Op 为 
“a7 人 ,4 (TT 
页 上 pr ( } (1) 
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计算 得 到 [sar 二 三 ”cpar 
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查 等 
| 分 界面 温度 为 Tp | 
图 2 两 材料 分 界面 温度 计算 流程 图 


Fig.2 Calculation flowchart of interface temperature between 


判断 


two materials 
Dh 
上 hn (7 oO 


因此 ， 车 一 一 > 则 Tp 的 取 值 正 
”Ar 


确 ， 否 则 调整 Tjp 的 取 值 重新 计算 ， 直到 得 到 正确 的 
分 界面 温度 。 同 理 可 求 得 N 型 电 偶 臂 分 界面 温度 Tv。 
P、N 电 偶 辟 的 平均 塞 贝 元 系 数 wp 和 ww 分 别 为 

人 wa)a7 + | ”aaa7 


Q, = Tp ZL (2) 
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其 中 xpj(D 、apz(D) 、aw(TD) 、axwz( 妃 分别 为 P_、N 
型 高 低温 材料 的 塞 贝 克 系 数 随 温度 变化 的 函数 。 


P 型 和 NN 型 电 偶 臂 的 总 内 阻 Rp、Rw 分 别 为 
[patDar + [ pn Da 图 
2 
PF pn Da + 六 md/ 和 


项 二 
A 
其 中 ppy()、pp2(D)、PN1(D)、pnz2(D) 分 别 为 P、N 
型 高 低温 材料 的 电阻 率 沿 电 偶 臂 方 向 变化 的 函数 ,4 
为 电 偶 辟 的 横 截 面积 

则 温差 电 单 偶 的 内 阻 Ri; 为 


Rin=RP+RN (6) 
温差 电 单 偶 产生 的 开路 电压 U 为 
RD 0) 
其 中 w 为 总 塞 贝克 系数 ， 由 a=ap-an 计 算得 到 。 
单 偶 两 端 连接 一 负载 R.， 则 电路 中 电流 厂 为 
ya a UV a al(7, 7) 8 
”二 名 


由 于 陶瓷 片 两 端 存在 温差 ， 冷 热 端 两 块 陶瓷 片 
在 工作 中 传递 的 热量 0 和 2. 分 别 为 


= (9) 
0. = (7 — TN, (10) 
其 中 
hd 
We (11) 
L.A. 
k= (12) 
后 和 .分 别 为 扫 端 和 准 端 向 交 片 的 六 导 ， i 
An、 ln 和 4.、 过 ss c 刀 分 别 为 热 端 和 冷 今 端 陶瓷 片 的 导 :; 热 


系数 、 横 截面 积 和 长 度 。 


根据 珀 尔 贴 定律 ， 当 电路 中 电流 为 时 ， 单 介 
热 端 吸 收 和 冷 端 释放 的 珀 尔 贴 热 9/、0; 分 别 为 
es (13) 
4 = olol, (14) 


民 据 焦耳 效应 ， 当 电路 中 电流 为 jy 时， 温差 电 
单 偶 产生 的 焦耳 热 2 为 


Wy) 一 Lo 人 ( 1 5) 


民 据 傅 里 叶 定律 ， 由 于 温差 电 单 偶 冷 热 端 温 差 
产生 的 传导 热 Ox 为 
0 = 而 (7 一 也 (16) 
其 中 
kf, 1 1 
Oo 0， Oo 0 (17) 


And NA 44 44 

Ko 为 温差 电 单 偶 的 热 导 ，4p 、4pPp 、4w 、4p 分 
别 为 P、N 型 高 温和 低温 型 材料 的 平均 导热 系数 ， 其 
数值 根据 各 段 材料 两 端 温度 和 材料 导热 系数 随 温度 
的 变化 求 得 。 
由 于 珀 尔 贴 效应 、 侍 里 叶 效 应 和 焦耳 效应 的 影 
响 根据 能 量 守 恒 原 理 ，O;、0QO, 可 表示 为 中 


WY es Wi 丰 Or A 0.58) 
= glo7, + 有 (人 — 7) — 0.570 8, (18) 
= 二 -了 刀 )1 


YQ, = WY, + Yr +0.50/ 


= aTT, + Kon, — 7) + 0.570 R, (19) 


一 (7. 人 Jr 
2 功率 和 效率 方程 


温差 电 单 偶 连 接 负载 及 后, 电路 中 的 电流 为 1， 
此 时 温差 电 单 偶 的 输出 功率 P 为 


P= 12 =[e -7)/(R, + 下 


1n 


20 

= a(n, -7 ME + RY 0 
温差 电 单 偶 的 热电 转换 效率 7 定义 为 

n=P/0, (21) 


将 式 (15)、(13) 代 入 式 (16)， 可 得 到 温差 电 单 偶 
的 工作 效率 为 


7 二 IR /Lalo7, 击 Al, | 7.)- 0. 5 (22) 


3 结果 与 分 析 


3.1 热 面 温度 的 影响 
在 该 分 析 中 ， 选 取 冷 面 温度 为 300K， 单 偶 长 度 


个 hinaV i 全 1 作 甘 日 车 
CNINAaAAIVD1T 项 | | 


设 为 3mm， 选 取 热 面 温 度 由 500K 逐渐 增 大 的 到 
800K， 计 算 分 段 温 差 电 单 偶 性 能 随 负载 的 变化 。 

图 3 为 不 同 热 面 温 度 下 温差 电 单 偶 性 能 随 负 载 
的 变化 。 随 着 负载 逐渐 增 大 ， 单 偶 的 输出 功率 和 转 
换 效 率 均 先 增 大 后 减 小 。 当 负载 电阻 与 温差 电 单 偶 
的 内 阻 相 等 时 ， 温 差 电 单 偶 的 输出 功率 最 大 ， 此 时 
电路 中 的 负载 电阻 定义 为 最 大 功率 电阻 。 温 差 电 单 
偶 达 到 最 大 转换 效率 时 的 负载 定义 为 最 大 效率 电 
阻 。 对 比 图 3a 和 图 3b， 可 以 发 现在 相同 条 件 下 最 大 
效率 电阻 略 大 于 最 大 功率 电阻 。 
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b Conversion efficiency 


图 3 不 同 热 面 温 度 下 温差 电 单 偶 性 能 随 负载 的 变化 


Fig.3 Variation of performance of the thermoelectric unicouple 


with load resistance for different hot surface temperatures 

随 着 热 面 温度 逐渐 增 大 ， 单 偶 的 最 大 输出 功率 
逐渐 升 高 ， 如 图 3a 所 示 。 这 是 由 于 随 着 热 面 温度 的 
升 高 ， 温 差 电 单 偶 热 端 温度 快速 升 高 ， 冷 热 端 温差 
逐渐 增 大 ， 单 偶 产 生 电 压 明 显 上 升 ， 而 内 阻 虽 增 大 
但 变化 幅度 较 小 ， 因 此 最 大 输出 功率 逐渐 升 高 。 

随 着 热 面 温度 的 升 高 ， 分 段 温 差 电 单 偶 的 最 大 
转换 效率 逐渐 增 大 ， 增 大 速度 逐渐 降低 ， 如 图 3b 所 
示 。 随 着 热 面 温度 的 升 高 ， 两 种 材料 的 分 界面 温度 
逐渐 升 高 ， 低 温 型 BizTes 材料 两 端的 平均 温度 逐渐 
高 于 BizTes 的 最 佳 效 率 温度 , BizTes 材料 的 效率 出 现 
负增长 ， 虽 然 高 温 型 材料 的 效率 仍 在 逐渐 增 大 ， 但 
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分 段 温 差 电 单 偶 的 效率 增 速 逐 渐 降低 。 


随 着 热 面 温度 的 升 高 ， 最 大 功率 电阻 和 最 大 效 
率 电阻 略微 增 大 。 随 热 面 温 度 的 升 高 ，Bi2Tes 的 电阻 
率 逐 渐 增 大 ,CoSbs 电阻 率 随 温度 变化 较 小 , 且 CoSb3 
的 电阻 率 明 显 高 于 BiyTe3, 从 而 使 分 段 温 差 电 单 偶 内 
阻 的 变化 介 于 两 种 材料 之 间 ， 更 接近 于 CoSbs 电阻 
率 的 变化 ， 而 最 大 功率 电阻 与 温差 电 单 偶 的 内 阻 相 
等 ， 因 此 最 大 功率 电阻 随 热 面 温度 升 高 略微 增 大 。 
3.2 冷 面 温度 的 影响 
在 本 分 析 中 ， 设 定 热 面 温度 为 650K， 冷 面 温度 
设 定 为 由 300K 逐渐 增 大 到 500K， 单 偶 的 长 度 设 为 
3mm， 计 算 不 同 冷 面 温度 下 单 偶 的 性 能 随 负 载 的 变 
化 。 
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图 4 不 同 冷 面 温度 下 温差 电 单 偶 性 能 随 电阻 的 变化 


Fig.4 Variation of performance of the thermoelectric unicouple 


with load resistance for different cold surface temperatures 
图 4 为 不 同 热 面 温度 下 温差 电 单 偶 性 能 随 人 负载 
电阻 变化 的 曲线 。 随 着 冷 面 温度 逐渐 增 大 ， 温 差 电 
单 偶 的 输出 功率 和 热电 转换 效率 均 降 低 ， 这 是 由 了 
随 着 冷 面 温度 的 升 高 ， 温 差 电 单 偶 冷 端 温度 快速 上 
升 ， 冷 热 端 温差 逐渐 降低 ， 因 此 最 大 输出 功率 逐渐 
降低 。 由 于 冷 热 端 温差 逐渐 降低 ， 热 端 吸 热量 也 逐 
渐 降低 ， 但 降低 的 速度 低 于 输出 功率 的 降低 速度 ， 
因此 最 大 转换 效率 也 逐渐 下 降 。 
3.3 温差 电 单 偶 长 度 的 影响 
在 该 分 析 中 ， 热 面 温度 设 定 为 650K， 冷 面 温度 
设 定 为 300K， 分 别 计算 分 段 温差 电 单 偶 在 长 度 为 
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3mm、2mm、lmm、0.5mm、0.25mm 时 的 性 能 随 负 
载 的 变化 。 

图 为 Sa 温差 电 单 偶 长 度 对 输出 功率 的 影响 ， 随 
着 温差 电 单 侦 长 度 由 3mm 逐渐 缩短 ， 分 段 温 差 电 单 
偶 的 输出 功率 均 逐 渐 增 大 ， 最 大 功率 电阻 逐渐 减 小 。 
这 是 由 于 随 着 温差 电 单 偶 长 度 减 小 ， 其 热 阻 减 小 ， 

冷 热 端 温差 减 小 ， 产 生 的 开路 电压 也 随 之 减 小 ， 同 
时 温差 电 单 偶 的 内 阻 也 减 小 ， 且 减 小 速度 大 于 开路 


电压 的 减 小 速度 ， 故 最 大 输出 功率 逐渐 增 大 。 
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b Conversion efficiency 
5 不 同 温差 电 单 偶 长 度 下 温差 电 单 偶 性 能 随 负 载 的 变化 


Fig.5 Variation of performance of thermoelectric unicouple 


with load resistance for different unicouple thicknesses 

5b 为 温差 电 单 偶 长 度 对 转换 效率 的 影响 ， 随 
着 温差 电 单 偶 长 度 由 3mm 逐渐 缩短 ， 温 差 电 单 偶 的 
最 大 转换 效率 逐渐 降低 。 这 是 由 于 随 着 温差 电 单 侦 
长 度 的 缩短 ， 输 出 功率 虽 逐 渐 增 大 ， 但 由 于 长 度 缩 
短 ， 冷 热 端 传 热 热 阻 大 幅 降低 ， 温 差 电 单 偶 热 端 吸 
热量 大 幅 升 高 ， 故 最 大 转换 效率 逐渐 降低 。 
3.4 热 表 面 对 流 传 热 系数 的 影响 

在 该 分 析 中 ， 冷 面 温 度 设 为 300K， 温 差 电 单 偶 
长 度 定 为 3mm。 热 表面 有 温度 为 800K 的 流体 流 过 ， 
计算 在 理想 条 件 〈 热 表面 对 流 换 热 系数 无 穷 大 ) 和 
热 表面 对 流 换 热 系数 为 1300W/AmK、1000WAnK、 
500W/m*KK、100W/mK 时 性 能 随 负载 的 变化 。 
图 6 为 热 表 面 对 流 传 热 系数 对 输出 功率 的 影响 ， 
随 着 热 表 面 对 流 传 热 系 数 由 100W/m KK 逐渐 增 大 ， 


分 段 温 差 电 单 偶 的 最 大 输出 功率 和 最 大 转换 效 


率 均 


逐渐 增 大 ， 最 大 功率 电阻 略微 增 大 。 这 是 由 于 


随 着 


热 表 面 对 流传 热 系 数 的 的 增 大 ， 热 面 温 度 快速 升 高 ， 


冷 热 端 温差 上 升 , 根据 3.1 的 分 析 可 知 ， 此 时 最 
出 功率 和 最 大 转换 效率 均 逐 渐 增 大 ， 最 大 功率 


大 给 
电阻 


也 随 之 增 大 。 同 时 随 着 热 表面 对 流传 热 系数 的 增 大 ， 


最 大 输出 功率 和 最 大 转换 效率 的 增 速 均 逐 渐 降 
最 大 输出 功率 和 转换 效率 逐渐 接近 理想 条 件 。 
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鲁 对 流传 热 系数 下 温差 电 单 介 
变化 


图 6 不 同 热 表 性 能 随 负 


由 出 


氏 ， 


载 的 


Fig.6 Variation ofperformance with load resistance for 


different convection heat transfer coefficients of hot surface 


4 结论 


本 文 对 分 段 温差 
析 ， 计 算 了 不 同 条 但 
换 效率 随 负载 的 变化 趋势 ， 主 要 结论 总 结 如 下 : 


TT 


分 段 温 差 电 单 偶 的 最 大 输出 功率 。 


BE 单 偶 的 性 能 进行 了 模拟 分 
下 温差 电 单 偶 的 输出 功率 和 转 


1) 提高 热 面 温 度 、 降 低 冷 面 温度 、 缩 短 温差 电 
单 侦 的 长 度 和 提高 热 表面 对 流传 热 系 数 均 可 以 提高 


单 偶 长 度 、 提 高 热 表 面 对 流 传 热 系 数 均 可 以 提高 分 


段 温差 电 单 偶 的 最 大 转化 效率 


3) 热 表 面 的 换 热效率 对 温差 电 单 偶 的 输出 功率 


和 转换 效率 有 很 大 影响 ， 在 设计 温差 发 电器 时 应 尽 


可 能 提高 热 表 面 的 对 流传 热 系数 ， 同 时 设法 增 大 热 


端的 换 热 面积 。 
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